Journal of Biogeography

AimThe assembly of montane plant communities through time is underlain by historical and abiotic factors. However, the extent of evolutionary connectivity between ancient highland ecosystems and surrounding lowlands remains unclear. Here, we investigate the evolutionary connections between the campos rupestres, a hyperdiverse and fragmented montane vegetation complex in eastern South America, and lowland biomes surrounding it: savannas, rainforests, and seasonally dry tropical forests. LocationEastern South America. Time periodCenozoic. Major taxa studiedFlowering plants. MethodsUsing phylogenetic beta diversity analyses for 13 angiosperm clades, we assess the degree of lineage dissimilarity between campos rupestres subregions and adjacent biomes. We also apply generalized dissimilarity modeling to determine the role of climate, soil, and geographic distance in shaping spatial patterns of phylogenetic composition. ResultsOur results reveal high lineage permeability between campos rupestres and surrounding biomes, with lineage sharing largely reflecting biome adjacency. This pattern is mainly driven by shared climatic conditions, which are the strongest predictors of phylogenetic dissimilarity. Main conclusionsWe highlight the importance of lineage exchange between lowland and montane environments for the assembly of highland floras. By showing that lineage movements across biome boundaries have been common over time and spatial scales, our study challenges the idea that ancient Neotropical mountains are isolated sky-islands. Instead, we emphasize the dynamic nature of montane plant diversity and the pivotal role of climate in shaping evolutionary connections between highlands and lowlands.

AimThe Andean paramo is the most biodiverse high-mountain region on Earth and past glaciation dynamics during the Quaternary are greatly responsible for its plant diversification. Here, we aim at identifying potential climatic refugia since the Last Glacial Maximum (LGM) in the paramo, according to plant family, biogeographic origin, and life-form. LocationThe paramo region in the Northern Andes MethodsWe built species distribution models for 664 plant species to generate range maps under current and LGM conditions, using five General Circulation Models (GCMs). For each species and GCM, we identified potential (suitable) and potential active (likely still occupied) refugia where both current and LGM range maps overlap. We stacked and averaged the resulting refugia maps across species and GCMs to generate consensus maps for all species, plant families, biogeographic origins and life-forms. All maps were corrected for potential confounding effect due to species richness. ResultsWe found refugia to be chiefly located in the southern and central paramos of Ecuador and Peru, especially towards the paramo ecotone with lower-elevation forests. However, we found additional specific patterns according to plant family, biogeographic origin and life-form. For instance, endemics showed refugia concentrated in the northern paramos. Main conclusionsOur findings suggest that large and connected paramo areas, but also the transitional Amotape-Huancabamba zone with the Central Andes, are primordial areas for plant species refugia since the LGM. This study therefore enriches our understanding on paramo evolution and calls for future research on plant responses to future climate change.

Background and aimsThe way plants cope with biotic and abiotic selective pressures determines their success in the colonization of remote oceanic islands, which ultimately depends on the phylogenetic constrains and ecological response of the lineage. In this study we aim to evaluate the relative role of geographical and ecological forces in the origin and evolution of the Madeiran ivy (H. maderensis). MethodsTo determine the phylogenetic placement of H. maderensis within the western polyploid clade of Hedera (three species), we analysed 40 populations (92 individuals) using genotyping-by-sequencing and including H. helix as outgroup. Climatic niche differences among the four study species were evaluated using a database with 706 records representing the entire species ranges. To test species responses to climate, a set of 19 vegetative and regenerative functional traits were examined for 70 populations (335 individuals). Key resultsPhylogenomic results revealed a nested pattern with H. maderensis embedded within H. iberica. Gradual niche differentiation from the coldest and most continental populations of H. iberica to the warm and stable coastal population sister to H. maderensis parallels the geographical pattern observed in the phylogeny. Similarity in adaptive traits is observed for H. maderensis and H. iberica. The two species show leaves with higher SLA, lower LDMC and thickness and smaller fruits than those of H. hibernica. ConclusionsAcquisition of the Macaronesian climatic niche and the associated functional syndrome in mainland European ivies (small fruits, leaves high SLA, and low LDMD and thickness) was a key step in the colonization of Madeira by the H. iberica/H. maderensis lineage, which points to climatic pre-adaptation as a driver of island colonization (dispersal and establishment). Once in Madeira, speciation was driven by geographical isolation, while ecological processes are regarded as secondary forces with a putative impact in the lack of further in situ diversification.

AimCommunity phylogenetics provides important information about the evolutionary and ecological factors help structure regional species assemblages. Here, we analyze phylogenetic diversity (phylodiversity) patterns among fish species in 97 sub-drainages of the Amazon basin, to evaluate the roles of historical and contemporary processes in generating and maintaining the exceptional richness and endemism of Amazonian fish species assemblages. LocationAmazon River basin TaxonFreshwater fishes MethodsUsing a large comprehensive database of freshwater fish species distributions, and a well-sampled molecular phylogeny of ray-finned (actinopterygian) fishes, we develop of multivariate statistical model to correlate estimated historical and contemporary environmental parameters with sub-drainage phylodiversity patterns. The model employs three phylogenetic metrics: i.e.: phylogenetic diversity (PD) sensu stricto, mean pairwise phylogenetic distance (MPD) between species capturing phylodiversity variation at older evolutionary timescales), and mean nearest taxon distance (MNTD) capturing variation in phylodiversity at younger evolutionary timescales. ResultsThe model recovered significant effects of elevation gradients, contemporary climate, habitat fragmentation, water types, and past marine incursions on assemblage phylodiversity patterns. The model also found significantly negative relationships among the three phylogenetic metrics, and between these metrics and distance to mouth of the Amazon, representing a West-East longitudinal gradient. Main conclusionsOur study revealed a highly non-random spatial and environmental distribution of our three phylogenetic diversity metrics across the 97 sub-drainages of the Amazon basin. Beyond significant regional effects of several environmental and historical drivers, we also found a significant West-East gradient of increasing phylogenetic diversity and phylogenetic relatedness, both patterns suggesting deeper evolutionary divergences among taxa located to the east, and more diverse, more recent radiations in the western sub-drainages. We conclude that western Amazonia can be seen as an evolutionary "cradle" of biodiversity for freshwater fishes in the Amazon basin as a whole. Significance StatementThis manuscript reveals spatial patterns of freshwater fish phylogenetic diversity and relatedness and explains its major contemporary and historical drivers in the Amazon basin. Amazon basin contains the highest freshwater biodiversity on Earth, as so investigate phylogenetic dimension of diversity is extremally relevant from the perspective of understanding the information on the evolutionary processes that had shaped Amazonian contemporary fish assemblages.

Evolutionary radiations on oceanic archipelagos (ROAs) have long served as models for understanding evolutionary and ecological processes underlying species diversification. Yet, diversity patterns emerging from ROAs have received relatively little attention from biogeographers, even though characterizing the effect of key geo-environmental factors on island clades species distribution could be important for unraveling diversification dynamics. In this study, we conducted a comparative analysis using island-specific species richness values for approximately one hundred ROAs across major oceanic archipelagos (mostly Hawaii, Canary Islands, Galapagos and Fiji) and taxa (vascular plants, invertebrates and vertebrates). Our aim was to determine whether (1) ROA species richness patterns scale as a function of key geo-environmental factors including island area, geological age, environmental heterogeneity (elevation and topographic complexity) and inter-island isolation, and (2) whether the magnitude of the effects of these factors varies across archipelagos and taxa. Our results identified elevation as a key driver of ROA species richness patterns on islands, supporting existing theoretical and empirical work that highlighted the central role of environmental heterogeneity in driving diversification on oceanic islands. As importantly, we found that the influence of geo-environmental factors varies across archipelagos and taxa, suggesting that unique archipelagic dynamics and biological traits together shape diversification differently. Our findings emphasize the value of applying biogeographical modeling at the resolution of individual radiations to improve our understanding of evolutionary processes on oceanic archipelagos.

AimTo evaluate the potential role of the orogeny of the Eastern Cordillera (EC) of the Colombian Andes and the Merida Andes (MA) of Venezuela as drivers of vicariance between populations of 37 tetrapod lineages co-distributed on both flanks, through geological reconstruction and comparative phylogeographic analyses. LocationNorthwestern South America MethodsWe first reviewed and synthesized published geological data on the timing of uplift for the EC-MA. We then combined newly generated mitochondrial DNA sequence data with published datasets to create a comparative phylogeographic dataset for 37 independent tetrapod lineages. We reconstructed time-calibrated molecular phylogenies for each lineage under Bayesian inference to estimate divergence times between lineages located East and West of the Andes. We performed a comparative phylogeographic analysis of all lineages within each class of tetrapod using hierarchical approximate Bayesian computation (hABC) to test for synchronous vicariance across the EC-MA. To evaluate the potential role of life history in explaining variation in divergence times among lineages, we evaluated 13 general linear models (GLM) containing up to six variables each (maximum elevation, range size, body length, thermoregulation, type of dispersal, and taxonomic class). ResultsOur synthesis of geological evidence suggested that the EC-MA reached significant heights by 38-33 million years ago (Ma) along most of its length, and we reject the oft-cited date of 2-5 Ma. Based on mtDNA divergence from 37 lineages, however, the median estimated divergence time across the EC-MA was 3.26 Ma (SE = 2.84) in amphibians, 2.58 Ma (SE = 1.81) in birds, 2.99 Ma (SE = 4.68) in reptiles and 1.43 Ma (SE = 1.23) in mammals. Using Bayes Factors, the hypothesis for a single temporal divergence interval containing synchronous divergence events was supported for mammals and but not supported for amphibians, non-avian reptiles, or birds. Among the six life-history variables tested, only thermoregulation successfully explained variation in divergence times (minimum AICc, R2 0.10), with homeotherms showing more recent divergence relative to poikilotherms. Main conclusionsOur results reject the hypothesis of the rise Andean Cordillera as driver of vicariance of lowland population because divergence dates are too recent and too asynchronous. We discuss alternative explanations, including dispersal through mountain passes, and suggest that changes in the climatic conditions during the Pliocene and Pleistocene interacted with tetrapod physiology, promoting older divergences in amphibians and reptiles relative to mammals and birds on an already established orogen.

Climate is well known as the main driver of species distributions. In this study, I focused on Geyik Mountains and surrounding areas to illustrate how complex topography and climatic conditions have shaped the distribution patterns of species/communities and therefore the biodiversity of the Mediterranean Taurus Mountains, one of the most biologically diverse areas in the Mediterranean Basin biodiversity hotspot. Accordingly, I used an ecological niche modelling approach, which has been widely used in recent biogeographic studies. I chose Taurus ground squirrels (Spermophilus taurensis) and coniferous forests as the representatives of high- and low-altitude species/communities, respectively. The results simply illustrate how complex topography and temperature and precipitation gradients have had a substantial role in shaping the distribution patterns of species/communities and therefore the biodiversity of the Mediterranean Taurus Mountains.

The tropical and subtropical Andes have among the highest levels of biodiversity in the world. Understanding the forces that underlie speciation and diversification in the Andes is a major focus of research. Here we tested two hypotheses of species origins in the Andes: 1. Vicariance mediated by orogenesis or shifting habitat distribution. 2. Parapatric diversification along elevational environmental gradients. We also sought insights on the factors that impacted the phylogeography of co-distributed taxa, and the influences of divergent species ecology on population genetic structure. We used phylogeographic and coalescent analyses of nuclear and mitochondrial DNA sequence data to compare genetic diversity and evolutionary history of two frog species: Pleurodema borellii (Family: Leiuperidae, 130 individuals; 20 sites), and Hypsiboas riojanus (Family: Hyllidae, 258 individuals; 23 sites) across their shared range in northwestern Argentina. The two showed concordant phylogeographic structuring, and our analyses support the vicariance model over the elevational gradient model. However, Pleurodema borellii exhibited markedly deeper temporal divergence ([≥]4 Ma) than H. riojanus (1-2 Ma). The three main mtDNA lineages of P. borellii were nearly allopatric and diverged between 4-10 Ma. At similar spatial scales, differentiation was less in the putatively more habitat-specialized H. riojanus than in the more generalist P. borellii. Similar allopatric distributions of major lineages for both species implies common causes of historical range fragmentation and vicariance. However, different divergence times among clades presumably reflect different demographic histories, permeability of different historical barriers at different times, and/or difference in life history attributes and sensitivities to historical environmental change. Our research enriches our understanding of the phylogeography of the Andes in northwestern Argentina.

AimThe effects of equilibrium and nonequilibrium processes are generally investigated using species richness on a single biological group. However, little is known about how these two classes of processes also affect trait diversity, considering multiple taxa within the same geographical template. Here, we evaluated which variables representing equilibrium (topography, climate, and primary productivity) and nonequilibrium (diversification rate and evolutionary time) processes best explain species richness and trait diversity of four clades of vertebrates within the same global hotspot. We also investigated how trait disparity has accumulated over time and whether there are congruent spatial patterns between groups. LocationAtlantic Rainforest Time periodContemporary. Major taxa studiedTerrestrial vertebrates. MethodsWe tested whether the spatial pattern of Functional Dispersion (FDis), richness, diversification rate, and evolutionary time of each group are correlated. We used a spatially explicit structural equation model to test how species richness and functional dispersion are influenced by variables representing equilibrium and nonequilibrium mechanisms. Additionally, we explored how trait disparity accumulated over time in the four groups. ResultsWe found that non-equilibrium proxies, evolutionary time and diversification rate, played a primary role in driving species richness and trait diversity, with elevation and climate variables having only an indirect effect species and trait diversity via diversification rate and evolutionary time. We found a congruent pattern of species richness among all groups, except among ectotherms. In contrast, the spatial distribution of evolutionary time was distinct for each group. Main conclusionsDespite nonequilibrium processes were more important for generating large-scale diversity patterns within the same geographical template, the interplay between evolutionary time and dispersal ability have disparately determined the assembly of communities.

AimBiodiversity is distributed unevenly among lineages and regions, and understanding the processes generating these global patterns is a central goal in evolutionary research, particularly in light of the current biodiversity crisis. Here, we integrate phylogenetic relatedness with species diversity patterns in four major clades of living tetrapods (amphibians, squamates, birds, and mammals) to approach this challenge. LocationGlobal. Time period300 million years ago - Present. Major taxa studiedTetrapods. MethodsWe studied geographic patterns of richness-corrected phylogenetic diversity (residual PD), identifying regions where species are phylogenetically more closely or distantly related than expected by richness. We explored the effect of different factors in residual PD: recent speciation rates, temporal trends of lineage accumulation, and environmental variables. Specifically, we searched for evolutionary and ecological differences between regions of high and low residual PD. ResultsOur results show heterogeneous spatial patterns of diversity dynamics across tetrapods. They reveal an overall negative relationship between recent speciation rates and residual PD, underscoring the role of recent speciation events in structuring current biogeographic patterns. Furthermore, we found differences between endothermic and ectothermic tetrapods in response to temperature and precipitation, highlighting the pivotal role of thermal physiology in shaping diversity dynamics. Main conclusionsGeographic patterns of diversity dynamics are heterogeneous across tetrapod clades and help us disentangle the evolutionary and ecological processes underlying them. By illuminating the multifaceted factors underpinning global diversity patterns, our study represents a significant advancement towards better understanding of how the present-day diversity of tetrapods was formed and how speciation rates influenced their species and phylogenetic diversity across clades and regions.

Under climate change, the status of invasive species is becoming complex, as taxa considered invasive in one region may become threatened in another, experiencing simultaneous expansion and decline across their range. In this context, mountain systems may act both as climatic refuges and/or as out-of-reach areas for these species. The Pyrenees, a longitudinal barrier separating the Iberian Peninsula from the rest of Europe, are increasingly exposed to alien plant species spreading from lowlands. Several invasive species have already reached high elevations in the massif, raising concerns about how many more may follow and to what extent under future climate change. This study assesses future invasion risks for 35 non-native plant species recognized as invasive on both sides of the Pyrenees. We project future dynamics of their bioclimatic niche suitability to identify which are likely to become widespread and where by 2100. Combining occurrence data with 19 bioclimatic variables at 1x1 km resolution, we model current bioclimatic niches and project its future distribution using an ensemble framework integrating five algorithms under four climate scenarios covering the 2021-2100 period. Projections indicate that 80 % of invasive plant species in the Pyrenees would not be climatic winners under future warming. The centroid of suitable climatic conditions is projected to shift upslope for more than 87 % of species. This upward movement would tend to stall around 2,000 m, as suggested by future hotspot patterns, effectively squeezing suitable areas and constraining future invasions within mountain zones, as climatic pressure increases at lower elevations. Such redistribution would not necessarily imply negative outcomes, and overlapping with endemic species would remain very limited. These results raise a broader question: could todays invasive species become tomorrows mountain flora?

Ecologists often agree on the importance of macroevolution for niche-mediated distribution of biological diversity along environmental gradients. Yet, macroevolutionary diversification and dispersal in time and space generate uneven geographic distribution of phylogenetic pools, which affects the imprint let by macroevolution on local species pools. In this article we introduce an individual-based simulation approach coupled to Approximate Bayesian Computation (ABC) that allows to parameterize the adaptation rate of species niche positions along the evolution of a monophyletic lineage, and the intensity of dispersal limitation, associated with the distribution of biological diversity between assemblages potentially connected by dispersal (metacommunity). The analytical tool was implemented in an R package called mcfly. We evaluated the statistical performance of the analytical framework using simulated datasets, which confirmed the suitability of the analysis to estimate adaptation rate and dispersal limitation parameters. Further, we evaluated the role played by niche evolution and dispersal limitation on species diversity distribution of Phyllostomidae bats across the Neotropics. The framework proposed here shed light on the links between niche evolution, dispersal limitation and the distribution of biological diversity, and thereby improved our understanding of evolutionary imprints on ecological patterns. Perhaps more importantly, it offers new possibilities for solving the eco-evolutionary puzzle.

Sky islands, mountain-top forests isolated by surrounding lowlands, offer unique opportunities to test how past and present landscapes shape species distributions. We examined the distribution of the Arizona gray squirrel (Sciurus arizonensis) across the Madrean Archipelago to test the constraint-based dynamic island biogeography (C-DIB) model, which posits that current occupancy in the sky islands reflects historical habitat size and connectivity. Using verified specimen records, we modeled climatically suitable habitats across four time periods: the Last Glacial Maximum (LGM), Mid-Holocene (MH), Present, and Future. For each mountain, we quantified suitable habitat area and estimated least-cost dispersal distances to assess both persistence and colonization potential. Our results suggest that species presence is best explained by LGM habitat metrics, which marginally outperformed models based on current conditions. Mountains that were large or well-connected during the LGM continue to support S. arizonensis, whereas historically isolated ranges remain unoccupied despite suitable contemporary habitat. These findings indicate a legacy of Pleistocene connectivity and reveal patterns of distributional disequilibrium. Furthermore, climatically suitable habitat for S. arizonensis has shifted both elevationally and geographically through time, reflecting long-term responses to climatic change. Together, these results emphasize the importance of protecting historically connected refugia, restoring riparian corridors that facilitate dispersal, and developing mountain range-specific management strategies that account for elevational shifts and potential downslope habitat recovery under future climate scenarios.

AimTo assess if and how species range size relates to range structure, if the observed geographic range properties can be retrieved from predicted maps based on species distribution modeling, and whether range properties are predictable from biogeophysical factors. LocationEurope Time periodCurrent Major taxa studied813 vascular plant species endemic to Europe MethodsWe quantified the size and spatial structure of species geographic ranges and compared ranges currently occupied with those predicted by species distribution models (SDMs). SDMs were constructed using complete occurrence data from the Atlas Florae Europaeae and climatic, soil and topographic predictors. We used landscape metrics to characterize range size, range division and patch shape structure, and analysed the phylogenetic, geographic and ecological drivers of species range size and structure using phylogenetic generalized least squares (pGLS). ResultsRange structure metrics were mostly decoupled from species range size. We found large differences in range metrics between observed and predicted ranges, in particular for species with intermediate observed range size and occupied area, and species with low and high observed patch size distribution, geographic range filling, patch shape complexity and geographic range fractality. Elevation heterogeneity, proximity to continental coasts, Southerly or Easterly geographic range positions and narrow ecological niche breadth constrained species observed range size and range structure to different extents. The strength and direction of the relationships differed between observed and predicted ranges. Main conclusionsSeveral range structure metrics, in addition to range size, are needed to adequately describe and understand species ranges. Species range structure can be well explained by geophysical factors and species niche width, albeit not consistently for observed and predicted ranges. As range structure can have important ecological and evolutionary consequences, we highlight the need to develop better predictive models of range structure than provided by current SDMs, and we identify the kinds of species for which this is most necessary.

AimTo investigate how past climate change has shaped the genetic diversity and demographic responses of tarantulas across latitudes, and to test whether climate-demography relationships vary with latitude. LocationGlobal; spanning tropical to temperate regions. TaxonTarantulas (family Theraphosidae). MethodsWe compiled mitochondrial Cytochrome oxidase I (COI) sequences for 48 tarantula species worldwide, including newly generated sequences, to estimate nucleotide diversity ({pi}) and Tajimas D. Species distribution models (SDMs) were constructed under present-day and Last Glacial Maximum (LGM) climatic conditions to quantify changes in habitat suitability since the LGM. Using generalized linear models (GLMs), we tested whether genetic and demographic metrics were associated with latitude and climate-driven habitat change, and whether their relationship varied with latitude. ResultsPairwise correlations among latitude, habitat change, and genetic metrics showed no significant associations. However, GLMs revealed a significant interaction: the effect of habitat suitability change on Tajimas D was strongly positive at high latitudes but negative or negligible at low latitudes. This indicates that demographic responses to past climate change varied latitudinally. Several high-latitude species showed genetic signatures of demographic expansion and range increase since the LGM. Main conclusionsOur results support the hypothesis that species at higher latitudes experience stronger demographic fluctuations due to historical climate change, aligning with Darwins early predictions. Moreover, patterns of demographic growth in temperate taxa suggest that some species may benefit from recent warming, consistent with Janzens climatic variability hypothesis. These findings demonstrate that climate-driven genetic and demographic responses in tarantulas are shaped by latitude, highlighting the importance of integrating phylogeography with ecological niche modeling to understand species resilience under climate change.

Speciation is a dynamic process shaped by the interaction of gene flow, geographic isolation, and ecological divergence. The Oncocyclus irises of the Southern Levant represent a young radiation of narrowly endemic species considered to be in the course of speciation. In this study, we used RAD-sequencing and single nucleotide polymorphism analysis across nine described species to investigate patterns of genomic divergence, gene flow, and local adaptation. Phylogenomic analyses revealed a mix of well-supported clades for some species, previously defined by their morphology and distribution, and non-monophyletic lineages, with several species exhibiting shallow divergence and shared genetic ancestry. We found evidence for gene flow and historical introgression between Iris petrana, I. atrofusca and I. mariae, while other cases of non-monophyly appear driven by incomplete lineage sorting. Both geographic distance (IBD) and environmental factors (IBE), mainly altitude, temperature, and aridity, were significantly associated with genetic structure, suggesting that local adaptation contributed to divergence following range expansion. Based on our findings we propose that the divergence of the Oncocyclus iris species in the Southern Levant supports a stepping-stone dispersion model, in which north-to-south dispersal was followed by local adaptation, and introgression in secondary contact zones. Overall, these findings highlight the complexity of speciation and the need for integrative approaches to study the interplay between historical divergence, contemporary gene flow, and environmental differentiation in shaping genomic patterns.

AimIdentifying the mechanisms driving divergence in marine organisms is challenging as opportunities for allopatric isolation are less conspicuous than in terrestrial ecosystems. Here, we aim to estimate a dated phylogeny of the squid family Loliginidae, and perform ecological niche analyses to explore biogeographic and evolutionary patterns giving rise to extant lineages in this group, with particular focus on cryptic species with population structure along the western Atlantic coast. LocationWorld-wide. TaxonClass Cephalopoda, Family Loliginidae MethodsWe used three loci to infer gene trees and perform species delimitation analysis to detect putative cryptic speciation events. We then estimated a dated species tree under the Bayesian multispecies coalescent and used it to reconstruct ancestral distributions based on the currently known ranges of the species. Also, we tested the hypothesis of niche divergence in three recently diverged species subpopulations of the northwestern and southwestern Atlantic Ocean by ecological niche modeling and niche overlap measurement from occurrence data. ResultsThe phylogenetic analyses confirmed the monophyly for the current twenty-six species of the Loliginidae family. Our ancestral area reconstruction and divergence estimation revealed the origin and geographical dispersal of loliginid lineages. Additionally, the phylogenetic analysis and the species delimitation analysis supported geographic structure within D. pleii, D. pealeii and L. brevis. The ecological niche models revealed unsuitable habitat in the immediately adjacent area of the Amazonian Orinoco Plume, yet suitable habitat characteristics beyond this area. Main conclusionsOur study allowed us to confirm the monophyly of all currently recognized species within the Loliginidae family and we corroborate the biogeographical origin being the Indo-Pacific region in the Cretaceous. We found a possible new cryptic lineage and show evidence of the Amazon-Orinoco Plume as an ecological barrier, which influenced the diversification of this particular group of marine organisms.

A long-standing question in ecology is how are species population distributed across space. The highest abundance has been hypothesized to be in the spatial or niche center, though mixed patterns from empirical studies has triggered a recent debate. Here we propose a conceptual framework based on environmental suitability and dispersal to interpret the mixed evidence. We demonstrate that the highest abundance could occur in the spatial center, in the niche center, or somewhere in-between the two centers, depending on the environmental setup and dispersal ability. We found that spatial and niche centers rarely overlap, suggesting the counteracting effect between the two factors, rather than reinforcement, is the norm in determining abundance patterns. The varied locations of highest abundance mirror the mixed evidence in literature, suggesting the "abundant-centre" and "abundant-niche centre" hypotheses are not mutually exclusive. This highlights the importance in understanding the biogeographic patterns through the lens of underlying mechanisms.

AimTo examine the support of two ecological diversity theories- The Ecological Limits Hypothesis (ELH) and the Niche Conservatism Hypothesis (NCH) - in explaining patterns of global shark diversity. LocationGlobal scale and two ecological realms: the Tropical Atlantic and the Central Indo-Pacific. Time PeriodPast 100 years Major Taxa StudiedWe examined 534 species of sharks and chimaeras, and we performed two subclade analyses on 272 species of ground sharks and 15 species of mackerel sharks. MethodsWe compared the species richness, mean root distance (MRD), and tree imbalance patterns to those simulated under the ELH and NCH with temperate and tropical centers of origin. We used sea temperature as a proxy for energy availability. We examined the importance of biogeographic history by comparing the model fits between two taxonomic groups, ground and mackerel sharks, and two geographic regions, the Tropical Atlantic realm and Central Indo-Pacific realm. ResultsThe ELH, temperate-origin model had the best fit to the global dataset and the sub-analyses on ground sharks, mackerel sharks, and the Tropical Atlantic. The NCH temperate-origin model provided the best fit for the Central Indo-Pacific. The {beta} metric of tree symmetry showed the best potential for differentiating between the ELH and NCH models, and the correlation coefficient for temperature vs MRD performed the best at differentiating between temperate and tropical origin of ancestors. Main ConclusionsThe global and subclade analyses indicate the ELH provides the best explanation for global scale shark diversity gradients even in clades with varying ecology. However, at the realm scale, biogeographic history has an impact on richness patterns. Comparing multiple metrics in relation to a simulation model provides a more rigorous comparison of these models than simple regression fits.

AimsSpecies distributional responses to climate change can depend on ecology and phylogeny. The degree of habitat specialism is potentially important because habitat generalists with wider distributions are assumed to be less sensitive to environmental changes compared to narrowly-distributed habitat specialists. Additionally, predicting which aspects of the changing climate impacts species with different habitat associations may be particularly challenging in complex environments such as those of tropical mountains. We explore the effects of climate change on potential distributions of Muscicapidae Flycatchers in parts of the Eastern Himalayan Mountains and the Indo-Burman range, and relate predictions to the level of ecological habitat specialism, while accounting for phylogenetic relatedness. LocationBhutan and Northeast India MethodsWe used Maxent to develop ecological niche models of 60 Muscicapidae Flycatchers under different climatic scenarios by collating presences from Global Biodiversity Information Facility, Bioclimatic variables from WorldClim, and elevation from ASTER. Species were scored as habitat specialists or generalists using the Species Specialism Index. Variables with high contributions to Maxent models were extracted to explore sensitivity to climate change based on habitat specialism while testing for phylogenetic signal using Phylogenetic Generalised Least Squares. ResultsMaxent models had the highest contributions from variable bio8 (mean temperature of wettest quarter) under present climate, and tmax (maximum temperature) under future climatic scenarios. Phylogenetic Generalised Least Squares revealed that habitat generalists had higher sensitivity to climate change than specialists. We did not detect strong phylogenetic signal in sensitivity to abiotic variables under all climate scenarios. Main conclusionsPotential distributions of Muscicapidae Flycatchers were sensitive to temperature variable in the month of the highest precipitation, and to maximum temperature. Potential distributions of habitat generalists were particularly sensitive to these abiotic variables, and those of habitat specialists less so. Sensitivity to abiotic variables did not show a pattern of phylogenetic niche conservatism. climate change, Muscicapidae Flycatchers, phylogenetic niche conservatism, habitat specialism, tropical mountains, Maxent